Gutzwiller Magnetic Phase Diagram of the Undoped t — t' — U Hubbard Model 
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We calculate the magnetic phase diagram of the half-filled t — t' — U Hubbard model as a function 
of t! and U, within the Gutzwiller approximation RPA (GA+RPA). As U increases, the system 
first crosses over to one of a wide variety of incommensurate phases, whose origin is clarified in 
terms of double nesting. We evaluate the stability regime of the incommensurate phases by allowing 
for symmetry-breaking with regard to the formation of spin spirals, and find a crossover to com- 
mensurate phases as U increases and a full gap opens. The results are compared with a variety of 
other recent calculations, and in general good agreement is found. For parameters appropriate to 
the cuprates, double occupancy should be only mildly suppressed in the absence of magnetic order, 
inconsistent with a strong coupling scenario. 

PACS numbers: 



I. INTRODUCTION 

An important issue in the Hubbard model is the na- 
ture of the metal-insulator transition, whether it is more 
Mott-like (driven by suppression of double occupancy, 
with no accompanying magnetic order) or more Slater- 
like (associated with magnetic order and a Stoner-factor 
instability). The question is rather subtle, since for in- 
stance a Mott phase can have a significant exchange cou- 
pling, which can lead to a parasitic magnetic order at 
low temperatures. Alternatively, critical fluctuations in 
a two-dimensional system will drive the magnetic order- 
ing temperature to zero (Mermin- Wagner theorem) while 
leaving behind a finite temperature pseudogap. Recently 
Tocchio, Becca, Parola, and Sorella (TBPS) 1 2 carried 
out a variational calculation of the t — t' — U Hubbard 
model at half filling, and found that the whole T = 
phase diagram, Fig. 1(a) is dominated either by paramag- 
netic phases or by phases with long-range magnetic order 
(solid green line), except for a small window of nonmag- 
netic insulator ('spin liquid' - dashed green line). Here we 
show that the ordered magnetic phase boundaries can be 
well reproduced by simpler Gutzwiller calculations, and 
that all of these instabilities are only weakly renormalized 
from the RPA values. These calculations are sufficiently 
simple that full allowance can be made for incommensu- 
rability [TBPS only studied antiferromagnetic order at 
(7r,7r) and (7r,0)], leading to a much richer phase dia- 
gram. The domain where TBPS found the spin liquid 
phase is characterized by a large number of competing 
phases, leading to potential frustration. 

Using a Gutzwiller approximation (GA), Brinkman 
and Rice (BR) 3 -' found a sharp metal-insulator transi- 
tion at a critical U = Ubr, where the effective mass 
diverges and the average double occupancy goes con- 
tinuously to zero. The BR line is Ubr = 8\Ek\, where Ek 
is the average kinetic energy per carrier below the Fermi 
energy Ep. While the sharp second order transition is 
now known to be an artifact of the simplified variational 
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FIG. 1: (Color online.) Phase diagram, as a function of U and 
t' , showing Ubr (red solid line), and calculations from Refi^ 
(green lines). Shown also are GA+RPA calculations Uga in 
which the momentum of the instability is restricted to be 
(7r,7r) [dashed blue lines] or (n, 0) [dot-dashed blue lines]. 



scheme^, Ubr signals a crossover to a regime of small rid, 
and hence can still serve as a measure of strong correla- 
tions. 



However, we will show in the present paper that the 
BR transition usually takes place at much larger U than 
a magnetic instability towards an incommensurate mag- 
netic state. We will give a detailed analyis of these insta- 
bilites in terms of double nesting and discuss the differ- 
ences between Hartree-Fock (HF) and the GA approach 
concerning the magnetic phase diagram. Before present- 
ing the corresponding results in Sec. Ill we briefly intro- 
duce the model and formalism in Sec. II and conclude 
our discussion in Sec. IV. 
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II. MODEL AND FORMALISM 

Starting point is the two-dimensional one-band Hub- 
bard model 



and Di denote the variational ('double occupancy') pa- 
rameters. 

We have also defined the spinor operators 



H 



t- r 



u 



(1) 



where Cj i<T (c| CT ) destroys (creates) an electron with spin 

a at site i, and m tCr = c ity Ci >a . U is the on-site Hub- 
bard repulsion and tij denotes the hopping parameter 
between sites i and j. In the present paper we re- 
strict to hopping between nearest (~ t) and next-nearest 
(~ t')ncighbors leading to a dispersion in momentum 
space e£ = — 2i[cos(A; x ) + cos(fc y )] — 4t' cos(k x ) cos(k y ). 

Our approach is based on a generalized GA^ supple- 
mented with Gaussian fluctuations (GA+RPApUn order 
to evaluate the magnetic instabilities. Since in the fol- 
lowing we will also calculate spiral states we use a spin- 
rotational invariant Gutzwiller energy functional as de- 
rived e.g. in Ref. [7] 



and the z-matrix 
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In the limit of a vanishing rotation angle ip the z-matrix 
becomes diagonal and the renormalization factors 
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reduce to those of the standard GA (= zq for a param- 
agnetic system). Spiral solutions are then computed by 
minimizing E GA with respect to a homogeneous rotation 
of spins with wave-vecor Q 



Sf = S* cos(QR J ) 
S? = S sin(QR i ). 



(3) 
(4) 



In order to compute the magnetic instabilites one can 
derive an equation similar to the Stoner criterion Ujj p = 
l/max{xo{q)} in HF+RPA. Here 



Xo(q) 



1 ^ n k+q.CT — ™k,CT 



N ^ 

k.<T 



k+q 



denotes the bare static susceptibility, and the maximum 
is taken over all g-values. 

In order to derive the corresponding condition within 
the GA one has to calculate the response of the system 
to an external perturbation which couples to the spin 
degrees of freedom. This can be achieved by expanding 
the energy functional Eq. [2] up to quadratic order in the 
(spin) density fluctuationa^which yields the generalized 
GA Stoner criterion 



max{U GA Xo(q)} = 1, 



(5) 



with an effective magnetic interaction 
U GA = [N q + M q [2E x + M q {El - E%) X o/4]] l z l (6) 

The parameters iV q and M q = z (z' — z'^ ) are defined 

in the Appendix of Ref [TT] and 



Ei = \.2 — —- 



k+q 



4) 



q : cr 



(7) 



k.CT 



As discussed in Refs. 1101111 the effective magnetic in- 
teraction Uga < U and saturates as the bare U — > oo 
whereas in HF+RPA, U can be arbitrarily large. As a 
consequence, for a given momentum q HF+RPA yields 
a magnetic transition at any doping whereas the corre- 
sponding instabilities within the GA are usually confined 
to a specific doping range. 

In the next section we compare our results with those 
of TBPS which are based on a Gutzwiller variational 
calculation (without making the GA). In their approach 
magnetic phases are based on HF groundstates while the 
spin liquid phase is based on a BCS ground state. The 
calculations go beyond simple Gutzwiller by including a 
Jastrow factor and a backflow correction. 
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FIG. 2: (Color online.) (a) Phase diagram obtained includ- 
ing the regime of incommensurate phases, comparing Uhf 
(brown line) and Uga (blue line). Shown also is the line of 
first order transitions to commensurate (n, tt) or (tt, 0) or- 
der (light green solid line) and the (n, n) to (n, 0) crossover 
line (light green dashed line). For ease in comparisons, the 
TBPS crossover line (green dotted line) from Fig. 1 is repro- 
duced. To illustrate the dominant g-vector (Fig. 2(b)), the 
HF+RPA points are color coded and numbered consecutively, 
although only some of the numbers are indicated, (b) Posi- 
tion of the dominant susceptibility peak in HF+RPA, color 
code and numbers match the values in (a). Note that some 
extra, metastable points are included, not shown in (a). Blue 
dotted line traces evolution of stable points, (c) na/rido as a 
function of t' , assuming a constant U = 8t (red line) or Wt 
(blue line). Here is the double occupancy, and n^o = 0.25 
is its uncorrelated value. Symbols indicate experimental dis- 
persion renormalization Zdisp as a function of t' , from Ref. 1141 
which should be an experimental measure ofnd/rido- Letters 
refer to L = Laa-zSr^CuCU; N = Nda-rzCe^CuCUia; B = 
BiaSraCaCuaOg. 



III. RESULTS 

Figure 1 compares the t' /t — U phase diagram of 
TBPS 1 ^ (green lines) with the Brinkman-Rice transi- 
tion U br (red solid line) and with the Gutzwiller transi- 
tion Uqa, restricted to only (tt, tt) [(tt, 0)] magnetic order 
(dashed [dotted] blue line). It is seen that while the spin 
liquid phase falls above Ubr, magnetic phases arise for 
much lower Us, and the commensurate GA transitions 
are in excellent agreement with Utbps- Since the phase 
boundaries depend only on \t'\, we illustrate only the sit- 
uation t' < appropriate to the cuprates. 

However, in general competing incommensurate phases 
become unstable first, due to Fermi surface nesting, 
Fig. 2(a), which compares the full Uqa (blue line) with 
Uhf (brown line). While the HF+RPA calculation over- 



estimates the stability of the magnetic phase, the over- 
estimate is not very largeP^ Moreover, in all cases the 
most unstable q vector (along these symmetry lines) is 
the same for the HF+RPA and the GA+RPA calcula- 
tions. Thus the main effect of the GA is to renormal- 
ize U — > Uqa < U, thereby reducing the range of the 
magnetic ordered phases. In Fig. 2(a) the HF+RPA cal- 
culations are coded with variously colored circles which 
match the points in Fig. 2(b), denoting the ordering q- 
vectors. These changes are associated with the evolution 
of the Fermi surface (FS) with t! , as discussed below. 

Since the experimental U's in cuprates fall in the range 
~ 6 — 8t, Fig. 1 suggests the cuprates are closer to 
Slater than to Mott physics. Within the BR model, the 
double occupancy at half filling is given by rid/rida — 
1 — (U/Ubr) 2 = m/m*, where — 0.25 is the un- 
correlated double occupancy and m* the renormalized 
mass. In Fig. 2(c) estimates for the double occupancy 
rid are plotted within the GA, for an assumed U = 8t, 
lOt, representative of the cuprates. The relatively small 
reduction of n d explains why the HF+RPA results are so 
accurate. Since ~ l/ m *j the present results are con- 
sistent with the small observ ed enh ancement of the effec- 
tive mass (circles in Fig. 2(c))p3EH. Thus, in cuprates the 
double occupancy n d is reduced not by Mott physics but 
by strong AFM fluctuations. In fact, the mean- field AFM 
suppresses double occupancy so well that Gutzwiller pro- 
jection on an AFM ground state actually increases double 
occupancyP 

In Fig. 2, it must be kept in mind that Uhf and Uga 
define the onset of the magnetic instability, via a Stoner 
criterion. As U increases beyond the onset, a finite gap 
is opened and the optimal q can change. At large U the 
entire Fermi surface is gapped, and the g's which produce 
the largest gaps are favored. 

In order to access this crossover to commensurate 
phases we compute the energies of spiral textures with in 
the spin rotational invariant extension of the GA^EiD as 
outlined in Sec. II. The most stable spiral is determined 
by searching for the spiral wave-vector Q m i„ which min- 
imizes the GA energy. For values slightly larger than 
Uga we consistently find the minimum of the energy 
landscape £ m in(Qmtn) close to the momenta of the in- 
stabilities as shown in Fig. [2^i. Upon further increasing 
U/t and depending on t' /t the momenta Q m i n shift ei- 
ther towards the diagonal (0,0) — > (tt,tt) or towards the 
line (tt, 0) — > (7T,7r). For a critical U/t one then finally 
finds a first order transition towards either Neel order 
Q = (tt, tt) or towards linear antiferromagnetic (LAF) 
order at Q = (tt, 0) [or equivalently Q = (0, 7r)]P This 
first order transition involves a topological transition of 
the Fermi surface, from a spiral phase with pockets to a 
fully gapped commensurate phase. 

In Fig. [2^ we show the first order line (light green) 
as an upper boundary for the incommensurate regime 
together with Uga as the lower transition (blue line) 
between incommensurate spin spirals and paramagnet. 
The boundary between Neel and LAF order (light green 
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dashed line) is close to the corresponding transition found 
by TBP9U2] (green dotted line). Note that the (charge- 
rotationally invariant) GA samples all possible incom- 
mensuabilities, so we confirm that the two ordered phases 
found by TBPS are indeed the only allowed ordered 
phases for large enough U. However, within the present 
scheme there is no energy gain for projected BCS wave- 
functions in the repulsive Hubbard model. Therefore, 
while our simplified scheme allows for a detailed deter- 
mination of magnetic phase boundaries we cannot access 
the spin liquid regime found by TBPS (cf. [2J;). We also 
would like to point out that besides spiral textures the in- 
commensurate regime may also contain spin density wave 
solutions with an associated small charge density modu- 
lation. However, for selected solutions we have checked 
that this has only a marginal influence on the first order 
transition line which in Fig. [2] is of the order of the line 
width. 

Figures 3 and 4 show how the FS evolves with t', and 
how this is reflected in the peak susceptibility. Fig. 3(a) 
displays a map of Xo(q) for t' — — 0.55£, showing that xo 
is dominated by a complex series of ridges, partly defining 
a diamond-shaped plateau centered at (tt, 7t). The origin 
of these structures is readily apparent from Fig. 3(b) us- 
ing the concept of nesting curves. For the generic case of 
two Fermi surface segments, a nesting curve can be de- 
fined as the locus of all points q = kpi — kF2> where kpi 
is a point on the ith FS, FSj, with the restriction that 
when FSi is shifted by q it is tangent to FS2. For the pa- 
rameters of Fig. 2(a,b) one has one Fermi surface and the 
nesting curve is simply given by q = 2kf- where k^? is the 
(anisotropic) Fermi wave vector. The case of two Fermi 
surfaces is discussed below. The dashed lines are exten- 
sions of the nesting curves folded back into the first BZ. 
It can be seen that all of the sharp structure in xo falls 
along this nesting curve, and the susceptibility peaks cor- 
respond to points where two branches of the curve cross. 
These values indicate q- vectors which nest the FS, and 
the crossing points correspond to double nesting, along 
two separate regions of the FS (see Fig. 4, below). 

Figures 3(a,b) illustrate the typical behavior at small 
\t'\ (appropriate for cuprates). From Fig. 2(a), the peak 
susceptibility is seen to be approximately commensurate 
at (tt, tt) for \t'\ < 0.23t [brown circles], then becomes ver- 
tically incommensurate (tt — S,tt) for 0.23 < \t'/t\ < 0.55 
[blue circles]. For 0.55 < \t'/t\ < 0.72, a diagonal in- 
commensurate (tt — 5, tt — 5) phase arises [green circles]. 
Competition between these two phases, corresponding to 
the points A and B in Fig. 3(a), is also found in hole- 
doped cuprate^l. The positions of A and B are readily 
determined, Fig. 3(b). Peak A lies along the zone bound- 
ary (q x a — tt) with 

. It — Ep 

q y a = 2arccos[ _ }, (8) 
while peak B follows the zone diagonal (q x = q y ) with 
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FIG. 3: (Color online.) (a) Plot of bare susceptibility xo in 
the first Brillouin zone, for t' = —0.55, at half filling, (b) Sim- 
ilar plot with nesting curve (solid line) and its folded replicas 
(dashed lines). (c,d) Similar plots for t' = —0.73, with addi- 
tional susceptibility peaks and nesting curves, as discussed in 
text. 



q x a = 2arccos 



(9) 



The only exception to this rule is the extended nearly- 
commensurate region near (tt, tt) for small \t'\. 

The topology of the band dispersion undergoes a dras- 
tic change at \t'\ = 0.5i. At that ratio the Van Hove sin- 
gularity (VHS) moves down to the bottom of the band, 
leading to a dispersionless band along the x and y axes. 
Beyond \t' /t\ = 0.5, the VHS moves away from the band 
bottom, but now the point T changes from the band min- 
imum to a local maximum, opening the possibility of a 
second FS section. At \t'\ = 0.71t, this section crosses 
the Fermi level at half filling, leading to a more com- 
plicated susceptibility map, Fig. 3(c,d), with a second 
q = 2ki?2 nesting curve (blue line) centered on T, as 
well as a q = k^2 — kpi nesting curve (green line). In 
all cases the nesting curves match the positions of sharp 
structure in xo- [For the inter-FS nesting, the q- value cor- 
responds to shifting one FS until it is tangent to the sec- 
ond.] This leads to a greatly increased number of nesting 
curve intersections, and the xo peak shifts, first, briefly, 
to a kpi-rmxed nesting curve intersection (orange line 
in Fig. 2(b)), then to a /cF2-mixed nesting curve (violet 
line). Figure 3 provides an example of double nesting, 
showing two competing (/-peaks, corresponding to A and 
B in Fig. 3(c). Here nesting vector A [green arrows] in- 
volves double nesting of the large FS [green FSs], while 
in B [orange arrows] one nesting involves the large FS 
[brown FS], but the other involves nesting between the 
large FS and the T-centered pocket [orange FS] . 
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FIG. 4: (Color online.) Fermi surfaces (red and blue lines) 
and the g-shifted versions, illustrating two (competing) exam- 
ples of double nesting. Data are for t' = — 0.73t, x — 0, as in 
Fig. 2 (c,d) [arrows]. 



(close to numerical results) and not qualitative (gain/loss 
of new phase) changes. Hence, the nesting curves in- 
troduced herein should provide valuable tools for deter- 
mining the optimal nesting vectors for many different 
correlated materials, and how these evolve with doping, 
temperature, impurity scattering, etc. 

In commenting on the calculation of TBPS, Becca, 
Tocchio, and SorellaP summarized earlier attempts to 
characterize the phase diagram of the model: "Remark- 
ably, all these numerical approaches give very different 
results for the ground-state properties of this simple cor- 
related model. In fact, there are huge discrepancies for 
determining the boundaries of various phases, but also 
for characterizing the most interesting non-magnetic in- 
sulator." Here, we have presented two additional phase 
diagrams, in HF and GA+RPA approximations. A key 
point is that with increasing U the first instability is 
generically to an incommensurate phase controlled by 
Fermi surface nesting. For larger U the FS is fully 
gapped, nesting is unimportant, and the only stable mag- 
netic phases are the two commensurate phases found 
by TBPS. For the rest our phase boundaries agree with 
TBPS except for the spin liquid phase which is beyond 
the capabilities of a simple mean field approximation. 



IV. CONCLUSIONS 



Increasingly, it is becoming clear that many of the com- 
plications of strongly correlated systems have to do with 
competing phases, whether leading to nanoscale phase 
separation, 'stripes', or frustration. However, part of 
the problem is a very incomplete understanding of phase 
competition even at weak coupling, and how the compet- 
ing phases evolve between weak and strong coupling. In 
this sense the magnetic instability of the two-dimensional 
Hubbard model provides an ideal case study, having a 
parameter space which is two-dimensional (q Xl q y ). In 
HF+RPA U is constant and the competition reduces to 
finding the maximum of xo (Stoner criterion). We find 
that the Gutzwiller correction leads only to quantitative 
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